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Surface magnetism of Rh(OOl) from LDA+U calculations 
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We report calculations indicating the presence of a surface magnetic moment for Rli(OOl), moti- 
vated by the detection of a finite moment by magnetic linear dichroism experiments. We show that, 
while the density functional with the local density or generalized gradient approximations (LDA and 
GGA) for exchange and correlation yields a non-magnetic ground state, the application of the GGA 
plus on-site Coulomb interaction U method predicts surface magnetism, thus offering a solution 
to the long-standing discrepancy between experiment and theory. The calculated moment on the 
outermost Rh atom increases with the strength of the effective on-site parameter UcH = U — J, for 
Ucf! > 1.2 eV, and is as large as 1.24 for UcS ~ 2.5 eV. 

PACS numbers: 75.70.Rf, 71.20.Bc, 71.15.Ap, 71.15.Mb 



The finding of a surface moment for a material which 
is non-magnetic in the bulk could have far-reaching con- 
sequences. It would give valuable directions for design 
of novel nanostructurcd materials and could have tech- 
nological implications for the magnetic storage industry 
in the future. So far, the only experimental evidence 
of a magnetic surface on a non-magnetic material has 
been obtained for the rhodium (001) surfaced Goldoni et 
alX measured the linear magnetic dichroism in 3o? core- 
level photoemission of Rh(OOl) surface. They applied an 
external magnetic field and reversed it by 180°, which 
caused a difference in the core level intensity and shape, 
indicating the presence of a magnetic moment. In a later 
study^S the same authors concluded that the rhodium 
(001) surface displays either a weak ferromagnetic order- 
ing or super-paramagnetism at 100 K. The conclusion is 
based on their finding that a small residual magnetic field 
is necessary for the dichroic signal, but, at the same time, 
the intensity of the effect docs not depend on the magni- 
tude of the applied residual field. Earlier spin-polarized 
photoemission study by Wu et al^ found a weakly ferro- 
magnetic Rh(OOl) surface at room temperature with an 
estimated moment per surface atom between 0.1 /is and 
0.2 

The experimental findings are in line with the expec- 
tation that the tendency towards magnetism increases 
near metal surfaces, because of the narrowing of the 
density of states which yields a Stoner enhancement in 
the susceptibility. The elements which are close to sat- 
isfy Stoner criterion in the bulk phase are likely to form 
magnetic ordering in a reduced dimensionality. In par- 
ticular, there are theoretical and experimental efforts 
directed towards investigation of surface magnetism in 
some Ad and hd transition metals. For a few decades now, 
many theoretical studies have examined the properties of 
rhodium surfaces^ and possibilities of surface magnetism 
in rhodium^iSiii^ Rhodium clusters have been investi- 
gated theorcticall}i2ii2iii and experimentallyii^ Overall, 
it has been found that the moment depends on the size 
and the geometrical structure of the cluster, with the 



calculated values of the moments somewhat larger than 
the measured value of 0.8 per atom for a ten-atom 
cluster^ Rh monatomic wires have also been shown 
theoretically to be magneticii^ Calculations of rhodium 
overlaycrs on Ag and Aui^*i^ yielded a moment on Rh 
atoms. Recently, ferromagnetic order has been found 
theoretically in bulk Rh at the optimum lattice constant 
in the body centered cubic lattice>i& So far, no realis- 
tic ab-initio study has found a ferromagnetic solution for 
the rhodium (001) surface. Cho and Scheffier— found a 
non-magnetic ground state which was practically degen- 
erate in energy with the ferromagnetic state, obtained as 
a solution of the constrained, fixed spin-moment calcu- 
lations, for moments below 0.5 Similar results have 
been obtained more recentlyi^ with an imposed moment 
limited to the surface layers. These calculations indicate 
a high susceptibility on the rhodium surface and sug- 
gest that surface structures with defects and reconstruc- 
tions might be magnetic. Cho and Schefffer— also solved 
a long-standing problem of anomalously large theoretical 
surface relaxations compared to the experimental data, 
by taking into account the thermal expansion of the sur- 
face, within the quasi-harmonic approximation. 

In this paper, we explore the possibility for mag- 
netism in rhodium (001) surface within the framework 
of the density-functional theory (DFT), using the gen- 
eralized gradient approximation (GGA) of the exchange- 
correlation potential, and the GGA plus on-site Coulomb 
interaction approach (denoted as LDA-f U from this point 
on). As the discrepancy between theoretical and ex- 
perimental results on the surface magnetism of rhodium 
persists for all previously attempted band-structure ap- 
proaches, it seems likely that the approximation used for 
exchange and correlation may be the reason for the non- 
magnetic solution. There are various examples in the 
literature, not limited to strongly correlated systems, for 
which the LDA+J7 was originally designed, where appli- 
cation of this method corrected the DFT solutions by 
yielding results in agreement with experiments^ The 
LDA+C/ is a method going beyond the LDA by treat- 
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ing exchange and correlation differently for a chosen set 
of states, in this case, the rhodium Ad orbitals. The se- 
lected orbitals are treated with an orbital dependent po- 
tential with an associated on-site Coulomb and exchange 
interactions, U and J. In the LDA, the electron-electron 
interactions have already been accounted for in a mean 
field way, and, therefore, one needs to apply a double 
counting correction. There are several existing versions 
of this correction>i2ii242& We apply the most commonly 
used one, introduced by Anisimov et ali^ which satisfies 
the LDA atomic-like limit to the total energy. It is known 
as the self-interaction corrected (SIC) LDA-I-J7 method. 

We solve the DFT equations using the WIEN2k 
implementation^ of the full potential linear augmented 
plane wave (FLAPW) method in a supercell geometry. 
We model the surface using an 11-layer slab with a vac- 
uum thickness corresponding to 6 layers. For the sur- 
face layer relaxation, we applied the previously calculated 
result of —1.4 %, compatible with the experimental find- 
ing of —1.16 ± 1.6 We used the optimized lattice 
constant of 3.84 A, sphere radius of 1.28 A, energy cutoff 
equal to 13.8 Ry and fc— point sampling with (22 x 22 x 
1) fc— points mesh in the full Brillouin zone (66 k-points 
in the reduced Brillouin zone). The calculations have 
been performed utilizing the GGA in the form given by 
Perdew, Burke and Ernzcrhofi^ From convergence tests 
with the number of fc— points and energy cutoff, we esti- 
mate the numerical accuracy of the energy difference be- 
tween magnetic and non-magnetic solution to be 2 meV 
and of the surface magnetic moment to the 0.02 fis- We 
have also assessed the effect of the topmost layer relax- 
ation on the energy difference and the surface moment. 
Changing the relaxation to -2.8 %, changes the surface 
magnetic moment by 0.05 fiB and the difference between 
magnetic and non-magnetic surface energy by 6 meV per 
surface atom within the LDA-I-J7. 

Considering the fact that the Stoner criterion for bulk 
rhodium is close to being fulfilled, we first checked for the 
effects of the inclusion, within the GGA, of the spin-orbit 
interaction in Rh(OOl) surface calculations, which was 
not included in previous calculations. It did not result 
in a surface moment. A similar conclusion was obtained 
by checking the influence of steps and line defects on the 
surface, within the GGA. The removal of every other line, 
or a pair of lines of surface atoms, parallel with the [110]- 
direction, gave only a non-magnetic solution. Structures 
with a few lines of surface atoms removed, simulating a 
step, also yielded no magnetic moment. 

In contrast, the LDA+J7 method did induce a mag- 
netic moment on the clean, defect-free surface. Ta- 
ble ^ gives the information on how the magnetic mo- 
ments in the surface and subsurface layers change as a 
function of the strength of the on-site parameter J/oS, 
given by: Ues = U — jMi^ We estimated the val- 
ues of the Coulomb and exchange parameters for bulk 
rhodium utilizing the Tight Binding Linear Muffin-tin 
orbital code2& in the atomic sphere approximation (TB- 
LMTO-ASA) and obtained U « 3.4 eV and J « 0.6 eV, 
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TABLE I: Magnetic moments on the surface (ms) and sub- 
surface (7714) layers, and difference of surface energies between 
magnetic and non- magnetic solutions, ASam-f, as a function 
of the effective Coulomb parameter UeS- 



i.e. UeS = 2.8 eV. Solovyev et alm^ calculated U and J 
using the same TB-LMTO code for rhodium impurities 
in Rb. They obtained a value oi U = 3.6 eV, J = 0.6 eV, 
t^cff = 3.0 eV for the monovalent Rh^+ impurity, which 
is not too different from the value we find. It is known 
that the calculated values of U and J tend to depend 
somewhat on the method used to calculate them and 
consequently have relatively large error bars (of the or- 
der of 1 eV). Therefore, the calculated values are to be 
used as a guidance towards a reasonable physical range 
of U values and an effort should be made to observe and 
analyze trends as a function of U. In Tabic H we re- 
ported the magnetic moments calculated for Ucs up to 
2.5 eV, together with the corresponding difference be- 
tween the magnetic and non-magnetic surface energies: 
Ai?surf = Ef^^f — E^ff . The non-magnetic energies were 
obtained using fixed-spin calculations with moment set 
to zero. 

For values of Ucs below 1.2 eV there is no ferro- 
magnetic solution. Starting with Ueff of 1.2 eV, the 
ferromagnetic solution has lower energy than the non- 
magnetic state. The surface energy difference for the case 
with Ucs = 2.5 eV is significantly larger than thermal- 
fluctuation energy at room temperature, thus indicating 
that the effect could be observable even at higher tem- 
peratures. With the increase of Ucs up to 2.5 eV, the 
magnetic moment on the surface and subsurface layers 
increases, and the ferromagnetic solution becomes more 
stable. We note that already for UeS ~ 1-5 eV, our cal- 
culated magnetic moments are one to two orders of mag- 
nitude greater than those obtained by Niklasson et alm^ 
from LDA+U calculations for Rh overlayers of compa- 
rable thickness on Ag(OOl), and do not disappear with 
increasing thickness^. 

The magnetic moments per atom are presented in 
Fig.^as a function of the layer position inside the slab for 
the case C/cff = 2 cV. The center of slab (denoted as layer 
0) has a vanishing magnetic moment and the values of the 
moment are clearly increasing towards the subsurface and 
surface layers (layers 4 and 5). We have also performed 
LDA-|-[/ calculations on bulk rhodium and found that, 
for all values of Ucs used in this work, the non-magnetic 
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solution is the stable state for the bulk. Considering the 
fact that the Rh Ad orbitals are not very localized, we 
performed a test using a different basis for the LDA + 
t/ in a pseudopotential plane-wave codei^ In this case, 
the localizcd-orbital basis consisted of the orthonormal 
atomic tight-binding orbitals, while in the FLAPW ap- 
proach, the orbitally dependent potential was applied to 
the atomic truncated orbitals (and was zero outside the 
muffin-tin radius). In the plane- wave scheme, we used 
an ultrasoft pseudopotential^ with a nonlinear core cor- 
rection and Perdew-Burke-Ernzerhof exchange and corre- 
lation potentiali^ The calculations gave similar results, 
with the values of magnetic moments somewhat larger 
than in the FLAPW code. For the sake of completeness, 
we did tests using the LDA for exchange and correlation 
and keeping all the other parameters of calculation the 
same. Just like in the case of GGA, the magnetic mo- 
ment on the surface atom increases and the ferromagnetic 
solution becomes more stable with increasing UcS, albeit 
with slightly smaller value of the moment on the surface 
Rh atomiS Finally, we also checked for the robustness 
of the Rh surface magnetism using another version of 
the LDA-f [/, with a different double counting correction, 
implemented in the WIEN2k code, the so-called mean- 
field correctioni^ For somewhat higher values of Ucs it 
also resulted in a ferromagnetic solutioni^ Hence, in all 
cases we obtain a stable surface magnetic solution, with 
a significant moment on the outermost atomic layer, for 
values of Ucs smaller than, or comparable to the cal- 
culated Ucs- Larger values of the on-site interaction Ucs 
induce thicker magnetic surface regions, requiring the use 
of thicker slabs. 
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FIG. 1; (Color online) Distribution of magnetic moments per 
atom inside a 11-layer Rh(OOl) slab. Layer is the center and 
layer 5 the surface of the slab. 
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FIG. 2: Spin density difference map (in e/A^), for the case 
UcH~2 eV, within a (100) atomic plane perpendicular to the 
surface. 



To gain a further insight into the microscopic arrange- 
ment of magnetic moments, in Fig. |21 we show the spin 
density difference map within a (100) atomic plane per- 
pendicular to the surface for the case Ucs = 2 eV. The 
spin density is very localized on the atoms of the surface 
and subsurface layers and is vanishing towards the cen- 
ter of the slab (bottom line of atoms in Fig.l^J- The spin 
polarized density on the outermost atoms derives mostly 
from dxy, d^z^-r^, dxz and dyz orbitals, with a contri- 
bution of ^ 30% from the d^y orbitals and of 21-22 % 
from each of the other three types of d orbitals. In fact, 
these four types of d orbitals dominate the local den- 
sity of states on the Rh surface atoms at energies around 
the Fermi energy (with comparable weights), both when 
UcS = and at finite UcS (2.5 eV). Inspection of Fig. [2 
also indicates that the states which contribute to the spin 
density map have predominantly antibonding character, 
which is consistent with the character of the d states 
found near the Fermi energy in bulk Rh. 

Intuitively, the magnetic ground state of the LDA-f J7 
method in this system can be understood by considering 
the effect of the on-site Coulomb interaction: it tends to 
favor the solution in which two spin states of the same d- 
orbital are separated in energy, with an energy separation 
roughly equal to U . To additionally check the reliability 
and stability of this method and its implementation, we 
applied it to (001) [or (0001) for hexagonal systems] sur- 
faces of some 4c? and 5o?-elements, namely Mo, Ru, W and 
Ir. All of these surfaces were found to be non-magnetic. 
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We used a relaxed geometry for all systems and, where 
appropriate, we took into account the surface reconstruc- 
tion (Mo and W) . In all cases we used the U ~ J values 
estimated in Ref. which should be an upper bound for 

f/eff. 

At present, a truly quantitative comparison between 
the experimental results and theory is not possible. Our 
solution, for which energy difference is greater than room 
temperature fluctuations (L'oft=2.5 eV), has a magnetic 
moment on the outermost Rh atom of 1.24 at K. 
Goldoni et alm^ could not give a measure of surface mo- 
ment, but concluded that their results arc consistent with 
weakly ferromagnetic or superparamagnetic surface. The 
only other experiment^ estimated the surface moment to 
be between 0.1 and 0.2 /i^. Experiments were done at 
100 K and room temperature, respectively. As we are un- 
able at present to include the temperature effects in our 
calculations, numerical comparison is not justified. It is 
also possible that a stronger ferromagnetic ordering in 
the experiment of Goldoni et alm^ would be detected, had 
the M states been probed directly. Their 3d core-level 
photoemission experiment probed the 4p polarization in- 
duced by the 4d moment. 

In summary, we have shown that a magnetic ground 
state for the rhodium (001) surface is obtained by us- 
ing the LDA+U approach. A stable solution with a 



significant magnetic moment on the outermost Rh atom 
is obtained for UcS smaller than, or comparable to the 
calculated C/cff, with a trend of increasing moment and 
stability of the magnetic solution with increasing UcS- 
We believe that our results on the rhodium (001) surface 
prove that the corrected non-local exchange and correla- 
tion in the LDA+J7 method is the key ingredient towards 
an accurate solution in density functional formalism and 
in agreement with experiment for this relatively weakly 
correlated system. It is our hope that similar approaches 
could be applied to other problems related to magnetism, 
even in the cases of relatively weakly correlated materi- 
als and including systems of lower dimensions^, in situa- 
tions where the DFT method itself cannot fully reproduce 
the experimental results. 
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